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ABSTRACT 

We present results from 4.8 GHz VLA and Global- VLBI observations of the northern half of 
the moderate FIR luminosity (median Lir = lO^^ '^^L©) COLA sample of star-forming galaxies. 
VLBI sources are detected in a high fraction (20/90) of the galaxies observed. The radio luminosi- 
ties of these cores (~10^^ WHz~^) are too large to be explained by radio supernovae or supernova 
remnants and we argue that they are instead powered by AGN. These sub-parsec scale radio 
cores are preferentially detected toward galaxies whose VLA maps show bright 100-500 parsec 
scale nuclear radio components. Since these latter structures tightly follow the FIR to radio- 
continuum correlation for star-formation we conclude that the AGN powered VLBI sources are 
associated with compact nuclear starburst environments. The implications for possible starburst- 
AGN connections are discussed. The detected VLBI sources have a relatively narrow range of 
radio luminosity consistent with models in which intense compact Eddington-limited starbursts 
regulate the gas supply onto a central super-massive black hole. The high incidence of AGN 
radio cores in compact starbursts suggests little or no delay between the starburst phase and the 
onset of AGN activity. 

Subject headings: galaxies: active — galaxies: starburst — infrared: galaxies — radio continuum: 
galaxies 
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1. Introduction 

It seems likely that super-massive black holes 
(SMBH) are ubiquitous in galactic nuclei. The 
nuclear stellar-velocity dispersio n and bulge-size 
to bl a ck hole mass correla ti ons dMagorrian et al 
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19981 : iMerloni et aP l2003t iHaring fc Ri>J \2004 



indicate that most galaxies with central bulges 
should contain SMBHs at their centers. The pres- 
ence of a SMBH is a necessary condition for AGN 
activity but only a small fraction of galaxies at 
any given time are strong AGN sources, imply- 
ing that a decisive factor determining AGN lumi- 
nosity is the gas feeding mechanism. While for 
Seyfert or LINER class AGN only relatively small 
black hole accretion rates are required (0.01 to 
0.1 M© year~^) gas originating in the galactic disk 
must lose angular momentum in order to be ac- 
creted onto the central black hole. Various mech- 
anisms to achieve such momentum transfer and 
m a.ss transport h ave been suggested (see review 
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napen 



2004 . 



There is increasing observatio nal evidence 



for a n AGN-starburst connection (iDavies et al 
20071: iHe ckman 200^ IPrieto et all l2005l l2007t 
Fathi et al.. ,2006.') . Plausibly such circumnuclear 



starbursts may play a role in removing angular 
momentum and so feedi ng the central black hole 
( Schartmann et al]l2010[ ). Alternatively both nu- 
clear starbursts and AGN could be independent 
consequences of a common gas transport mecha- 
nism into the centers of galaxies with no direct 
causal connection between the two phenomena. 
Understanding in detail the phenomenology of 
starburst and AGN activity, in particular whether 
one phase begins before the other, is important 
to test the nature of the starburst- AGN inter- 
relationship. 

Most detailed observational work so far on 
the starburst- AGN connection has studied star- 



burst activity in samples of AGN (jDavies et al 



20071 ). Here we adopt the alternative approach 



of searching for AGN activity in a far infrared- 
selected sample of star-forming galaxies known 
as COLA (Compact Obj ects in Low-power AGN, 
ICorbett etalll2002i [20031 see §2 for sample defini- 
tion). Our principal means of identifying AGN 
activity in this sample is by searching for high 
brightness temperature VLBI radi o cores. Classic 
optical emission line diagnostics (j Baldwin et al 



■Veilleux fc Osterbrockl Il987l : iKewlev et al 
Kauffmann et al.ll2003[ ) can be unreliable foi 



such searches as an AGN may be obscured by 
dust in star-forming galaxies. High brightness- 
temperature radio emission is unambiguous evi- 
dence of an AGN, however a non-detection does 
not exclude an AGN because of the existence 
of both radio lou d and radio quiet objects (see 
Zamfir et al.ll2008L for a recent review); it follows 
that VLBI searches can only give a lower limit 
to the AGN fraction. It should be noted that 
while radio supernovae (SNe) can also generate 
compact radio emission they have a maximum lu- 
minosity which dist inguishes them from AGN (see 
Corbett et al1l2002i and SJHU). 

We have already published radio observations 
l|Corbett et al.ll2002^) and optical spectroscopy re- 
sults Corbett et _aL (2003 ) for the S < 0° (COLA- 
S) part of the sample. Here we present radio and 
optical observations of the Northern half {6 > 0°) 
of the sample (COLA-N). This paper is organized 
as follows. In ij2]we give an overview of the COLA 
project and previous observations. In §3 we de- 
scribe our VLA, VLBI, optical spectroscopy and 
other observations and give an overview of the 
data reduction process. Our observational results 
are presented in §4 and our discussion in §5. Fi- 
nally in §6 we list our major conclusions. 

2. The COLA Project and Previous Ob- 
servations 

The COLA (Compact Objects in Low-power 
AGN) project has the primary goal of determining 
the relationship between AGN and other galactic 
properties (such as galactic structure, degree of 
interaction and star formation activity) in a mod- 
erate luminosity FIR selected sample. The sample 
contains all the galaxies in the IRAS Point Source 
Catalog (jBeichman et al.lll988r ) with (a) flux den- 
sities at 60yum greater than 4Jy, (b) heliocentric 
velocities between 3500 and 7000 km s^^ (to min- 
imize Malmquist bias) and (c) galactic latitude 
6|>10° (to avoid confusion with galactic objects). 
The resulting sample has 217 galaxies of which 110 
are located in the northern hemisphere (COLA-N) 
and 107 in the Southern hemisphere (COLA-S). 

The COLA sample is selected at the same IR 
wavelength as the Revised Bright Galaxy Sam- 
ple (RBGS see ISanders et al.l (|2003[ )) but with a 
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flux limit which is somewhat smaller (4 Jy versus 
5.24 Jy) plus an additional redshift range selec- 
tion. It follows that there is a significant over- 
lap between the two samples with 65/110 (59%) 
of the COLA-N sources also being members of the 
RBGS. 

Table 1 provides a summary of the COLA-N 
sample (for a similar table describing the south- 
ern sample see lCorbett et al. 2002). The luminos- 
ity distance Dl given column 7 for each source 
was calculate d in NED assuming the three attrac- 
tor model of iMould et aD (|2000[ ) and the cosmo- 
logica l parameters derived from 5 year WMAP 

= 0.276, 
Details on 



data (jHinshaw et al 



200i), 



-1 



0.726 and Ho = 70.5 kms 
the calculation of Infra-red and radio luminosi- 
ties presented in Table 1 are given in iJ3.3l Ta- 
ble 1 shows that most COLA sources have bolo- 
metric IR luminosities close to the boundary defin- 
ing Luminous Infra- Red galaxies (at Lui—IO^^Lq 
Sanders fc Mirabell ( 1996i) ) with a median IR bolo- 



metric luminosity for COLA-N of 10" °1Lq. All 
but one source has a bolometric IR luminosity in 
the range IQiO-S < Lir,/Lq <10" '^, the exception 
being the nearby ULIRG Arp 220. 

Low resolution 1.4, 2.5 and 4.8 GHz contin- 
uum observations of COLA-S sources have been 
made using the Australia Telescope Compact Ar- 
ray (ATCA). These radio observations were com- 
plemented with high resolution single baseline 
2.3 GHz snapshot observations obtained using the 
Australian Long Baseline Array (LB A). Optical 
spectroscopy for a large fraction of the southern 
sources was also obtained using the Dual Beam 
Spectrograph (DBS) on the 2.3 m telescope at the 
Mount Stromlo and Siding Springs Observatory. 
Detailed descriptions of the COLA- S radio and op 



tical observ ations can be found in ICorbett et al 
(|2002l |2003[ ) respectively. 



Corbett et aD (|2002l ) detected 9 out of the 105 
COLA-S galaxies in their LBA observations. Of 
these, 8 showed radio emission stronger than that 
predicted by the IR-radio correlation. Moreover 7 
of the 14 sources with a total radio flux density 
more than la above that predicted by the FIR- 
Radio correlation were LBA detections. Based on 



these results ICorbett et al.l (|2002l ) concluded that 
sources detected with the LBA exhibited a sta- 
tistically significant radio excess relative to the 
non-detections. This excess persisted even after 



the subtraction of the radio emission from the de- 
tected core implying the existence of associated 
AGN powered diffuse radio emission from struc- 
tures larger than ^15 pc. 



The COLA-S optical data (jCorbett et al.ll2003[ ) 
showed that in the southern sample the frac- 
tion of galaxies optically classified as Seyfert was 
■~15%. Of these only 55% were detected with the 
LBA suggesting the there may be two populations 
of Seyferts, one population with extended radio 
structures and compact radio cores and the other 
without. In this paper we revisit these tentative 
conclusions in light of the larger number of long- 
baseline detected sources in COLA-N. 

3. Observations 



3.1. VLA Observations 

Snapshot radio continuum observations for the 
COLA-N sources were obtained using the VLA at 
4.8 GHz in two observing epochs each of 24 hours 
in length (see Tabled . Sources with 0° < S < 30° 
were observed on 1998/06/23 in the VLA BnA 
configuration, while sources with 6 > 30° were ob- 
served on 2003/07/19 in the A configuration. In 
addition, to ensure that all compact sources were 
observed with the A-array, 12 sources found to 
be compact at 4.8 GHz from the first epoch were 
re-observed in the second epoch. The data were 
reduced using AIPS in a standard manner. Ra- 
dio maps of size 50"x50" centered at the IRAS 
positions were produced for all the sources using 
natural weighting. Typical FWHM beam-sizes at 
4.8 GHz were ~0.4" and 0.9"x0.4" for A and BnA 
configurations respectively. The typical rms noise 
levels were 75/iJy beam"-'^ 

3.2. VLBI Observations 

VLBI 4.8 GHz snapshot observations at data- 
rate 1 Gbit s~^ per station were conducted for 
the COLA North sample using the most sensitive 
telescopes available in the European VLBI Net- 
work (EVN). The first epoch on Feb 27/28th 2005 
used Eb (Effelsberg), Wb (Westerbork) and Ar 
(Arecibo) for all sample sources within the decli- 
nation range 10° < 5 < 40° (defined by the Ar 
zenith-distance restrictions and minimum eleva- 
tion at the European stations) with ~15 minutes of 
on-source time per target giving ct=25 /iJy on the 
Ar-Eb baseline. The second epoch on May 6/7th 
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2005 used Eb, Wb and Jbl (Jodrell Bank, Lovell 
telescope) for all the sources with 0° < i5 < 10° 
and S > 40° with ^15 minutes of on-source time 
per target giving noise a =56 ^Jy on the Eb-Jbl 
baseline. Gaps in the schedule were filled by re- 
observing 21 sources from the first epoch. In both 
epochs strong nearby phase reference and am- 
plitude calibrators (>500mJy) were observed for 
10 minutes every ~lhour. The Eb-Wb baseline 
(266 km) which is common to both epochs has a 
fringe spacing of 46.5 mas at 4.8 GHz which is sim- 
ilar to the 47.4 mas fringe-spacing of the ATCA- 
Tidbinbilla baseline (566 km) at 2.3 GHz used by 
Corbett et al. I (|2002h to observe the southern sam- 
ple. Because of scheduling restrictions and bad 
data we only obtained useful VLBI observations 
for 90 sources. 

AIPS was used for the initial data reduction 
and calibration stages with instrumental phase 
and delay being determined toward strong calibra- 
tors. Monitored system temperatures and gains 
were used to calibrate amplitudes. An adjustment 
of ~30% was required on Jbl to obtain consis- 
tent calibration on sources known to be unresolved 
from the VLBA calibrator database. Our final es- 
timate is that quoted baseline flux densities are 
accurate to ^-^15%. 

Fringe fitting was performed in AIPS on all 
the calibrators and then these solutions were ap- 
plied to the interleaved target source scans. Af- 
ter this we exported the data to our own software 
within which Delay-Rate maps were made for each 
baseline/target source. The size of these delay- 
rate maps were set such that a compact source 
would be detectable up to 20 " from the correla- 
tion position. Finally the Delay-Rate maps were 
transformed into RA-Dec offset maps relative to 
the phase center. These maps show only features 
whi ch are smalle r than the baseline fringe spacing 
(see lParrall2007l ). We adopted a detection thresh- 
old of 6 times the measured r.m.s. noise. The high 
sensitivity and small fringe spacing of the Eb~Ar 
baseline resulted in a very small uncertainty on 
the phase center offsets of the detections. These 
VLBI position offsets are shown as crosses on the 
corresponding VLA map in Figure [2] 

3.3. Archival Data Radio and IR data 

In addition to our own radio data our analy- 
sis made use of source 1.4 GHz total flux densi- 



ties t aken from the NVSS catalog (|Condon et al 
19981 ). Infra-red total flux densities were taken 



from the IRASPoint Source Catalog (jBeichman et al 
1988I ). In order to calculate source infrared 
to radio (1 . 4 GHz ) ratios {qia) as defined by 
Helou et"aLl ()l985[ ) we calculated for each source 
an equivalent FIR flux density from S'fir — 1.26 x 
10-^(2.58/60 + /ioo)/3.75 x lO^^ Wm^^Hz-i 
where /go and /loo are the cataloged 60 and 
100 /im /i?j4S'flux densities. For comparison with 
COLA-S and other data sets for each source a 
FIR spectral luminosity (see Table [U column 9) 
was calculated assuming isotropic emission. Fi- 
nally bolometric IR (8/xm -1000/xm) luminosities 
in solar units were calculated (see Table [U column 
7) from the IRAS flux densities usi ng the relation 
given bv ISanders fc Mirabell (|l996[) . 



3.4. Optical spectroscopy observations 

Optical spectroscopy observations were taken 
using the 60 inch CfA telescope at Mt Hopkins, 
Arizona between February and October 2001 us- 
ing the FAST spectrograph. A slit width of 3" 
was used, corresponding to 0.9 kpc at the median 
sample source distance. A grating of 300 lines/mm 
was used to disperse light between 4000 and 7000A 
and gave a spectral resolution of 6 A. A single expo- 
sure of between 900 and 1200 seconds was made for 
each target. Spectrophotometric and smooth stan- 
dards (often the same stars) were observed during 
each night and used to flux-calibrate and correct 
for atmospheric absorption features in the spectra. 
We were able to measure useful line ratios for 61% 
of the sample. 

Data reduction was carried out in the standard 
manner using IRAF. The frames were de-biased 
and flat-fielded. Galaxy spectra were extracted 
using an aperture centered on the position of peak 
flux across the slit. If more than one peak was seen 
in the spatial direction a spectrum was extracted 
at the position of each peak. Sky-subtraction was 
then carried out and wavelength calibration done 
using an Ne-Ar arc lamp. The observations were 
then corrected for atmospheric extinction and flux 
calibrated using the spectrum of a spectrophoto- 
metric standard. Finally atmospheric absorption 
features in the target spectra were removed using 
the smooth spectrum standard observed at a sim- 
ilar airmass to the target. 

The line fluxes were measured by fltting (mul- 
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Table 1 
COLA North Sample 



IRAS 


Other 




R.A. 




DccL 


VhoI 


Dl 


log LiR 


log Ll.4 


log Lfir 


91.4 


Name 


Name 




J2000 




J2000 


[km s"^] 


[Mpc] 


[i0] 


[W Hz"^] 


[W Hz"^] 




(1) 


(2) 




(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


00005+2140 


MRK334 


00 


03 09.74 


21 


57 36.46 


6579 


93.4 


11.07 


22.46 


24.74 


2.27 


00073+2538 


N23 


00 


09 53.58 


25 


55 26.40 


4566 


64.7 


11.09 


22.57 


24.82 


2.25 


00506+7248 




00 


54 04.00 


73 


05 11.70 


4706 


69.3 


11.49 


22.82 


25.21 


2.40 


00521+2858 




00 


54 50.27 


29 


14 47.60 


4629 


65.2 


10.88 


22.39 


24.62 


2.23 


00548+4331 


N317 


00 


57 40.58 


43 


47 32.30 


5429 


77.2 


11.16 


22.66 


24.92 


2.26 


00555+7614 




00 59 15 


76 30 52 


4739 


70.1 


10.91 


22.53 


24.66 


2.13 


01503+1227 




01 


52 59.52 


12 


42 27.90 


4558 


63.3 


10.91 


22.36 


24.65 


2.29 


01519+3640 




01 


54 53.93 


36 


55 04.40 


5621 


79.2 


11.02 


22.14 


24.74 


2.60 


01555+0250 


ARP126 


01 


58 05.26 


03 


05 00.90 


5431 


75.6 


10.94 


22.47 


24.68 


2.20 


01556+2507 




01 


58 30.63 


25 


21 36.90 


4916 


68.8 


10.99 


22.36 


24.75 


2.39 


01579+5015 




02 


01 09.65 


50 


30 25.50 


4875 


69.5 


10.85 


22.29 


24.55 


2.26 


02071+3857 


N828 


02 


10 09.52 


39 


11 25.30 


5374 


75.8 


11.33 


22.85 


25.07 


2.22 


02080+3725 


N834 


02 11 01 


37 39 58 


4593 


64.7 


10.95 


22.49 


24.69 


2.20 


02152+1418 


N877 


02 


17 58.62 


14 


32 25.90 


3913 


54.2 


10.98 


22.58 


24.74 


2.17 


02208+4744 




02 


24 08.00 


47 


58 10.68 


4679 


66.5 


11.09 


22.51 


24.84 


2.33 



Columns — (1): IRAS Name (from Point Source Calalog, excepting F07258+3357 which is taken from the Faint Source 
Catlog). (2): Other name. (3)-(4): Right Asccncion and Declination (from VLA observations, unless undetected in which 
case IRAS positions are given in italics). VLA positions arc conservatively estimated to be accurate to 0.3" and IRAS 
positions to 20". (5): Heliocentric velocity (cz) from NED. (6): Luminosity dist ance (see Section 2). (7): FIR. luminosity 
(see Section 3.3). (8): 1.4 GHz luminosity calculated from the NVSS catalog JCondon et al.||1998ll . (9): FIR spectra l 
luminosity (see Section 3,3). (10): qi, — log(LFiR/-£'i.4) • The full version of this table is given in Table|l (Full Version')] 



Table 2 

Summary of Radio Observations 



Date 


Array 


Frequency 


Beamsizc^ 


Sources 


Typical RMS 






[GHz] 


[mas] 




[^Jy beam"^]'^ 


98/06/23 


VLA BnA 


4.9 


900x400 


0° < (5 < 30° 


75 


03/07/19 


VLA A 


4.9 


400 


6 > 30° and CC ^ 


75 


05/02/27 


Eb-Wb-Ar'^ 


4.9 


1.8 


10° <S < 40° 


25 


05/05/06 


Eb-Wb-Jbl'^ 


4.9 


18.1 


0° < (5 < 10° and S > 40° 


56 



'^CC^Sources known to be compact at 4.9 GHz from the first VLA epoch. 

''For the VLA observations this is the typical FWHM of the restoring beam. For the VLBI 
observations this is the typical fringe spacing of the longest baseline. 

'^Eb=EfFelsberg, Wb=Westerbork, Ar=Arecibo and Jbl=Jodrell Bank (Lovell). Baselines 
lengths are Eb-Wb=266 km, Eb-Jbl=699 km and Eb-Ar=6901 km. 

'^r.m.s noise has units of /iJy beam^^ for VLA observations but for VLBI observations refers 
to noise on most sensitive baseline in /iJy. 
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tiple) Gaussians to the emission lines using 
NGAUSSFIT in IRAF. The continuum flux was 
subtracted using a linear fit to small portions of 
the continuum either side of the emission line. 
The Ha (A 6562A) and [Nil] doublet (AA 6583, 
6548A) were often overlapping, in these cases line 
profiles were de-blended by simultaneously fitting 
a separate Gaussian for each line. When Ha was 
present in absorption as well as emission both 
the emission and absorption features were fitted 
simultaneously with Gaussians. 

4. Results 

4.1. Sample overall Radio/FIR ratios 

In a FIR-selected sample of galaxies such as 
COLA, it is expected that in most sources star- 
formation will be the dominant mechanism pow- 
ering the FIR emission. For such sources radio 
and FIR luminosities and flux d ensities a re ex- 
pected to be tightly correlated (see Condonlll992.) . 
The total 1.4 GHz and FIR luminosities for our 
sample are plotted in Figure[T]-Left. Except for 
one outlier (3C84) which is unusually radio lumi- 
nous, sources tightly follow the expected correla- 
tion. We calculated for each source the logarith- 
mic ratio o f flux densities q^A = log(S'FiR/S'i.4) as 
defined bv IHcIou et all (|l985[ ). Excluding 3C84, 
the mean value of 91.4 is 2.34 ± 0.01 (and the me- 
dian 2.31 ± 0.01) which is in good agre ement with 
the m ean 91.4 — 2.34±0.01 measured bv lYun et al 



((2001') from a sample of 1809 galaxies (solid line in 
Fi gure[I]). In contrast t he median value for COLA- 
S (|Corbett et alll2002( ) is gi.4 = 2.43 ± 0.03 which 
is significantly (4cr) larger than the COLA-N me- 
dian value. 

To try to resolve the above discrepancy we 
compared NVSS and COLA fluxes for the 56 
COLA South sources also contained in the NVSS 
catalog (i.e. those with S > —40°). We found 
NVSS flux densities for these sourc es which were 
12% l arger than those measured by ICorbett et al. 
(I2OO2I ). resulting in qia = 2.36 ± 0.05 for ATCA 



fluxes and 2. 31 ib 0.04 for NVSS fluxes. It has 
been n oted by Prandoni et al. J 2000h. Norris et al 



(|2006l ) and iKellermann et al.l (|2008h that fo: 



sources with flux densities < 15mJy NVSS flux 
densities tend to be ^ 10% larger than ATCA 
measurements. Although COLA total flux densi- 
ties are about a factor of two above this flux limit 



the effect is neither understood nor well-defined 
and so this scale difference may contribute to the 
apparent difference of median gi.4. We conclude 
that the difference in gi.4 between COLA-N and 
COLA-S is likely caused by (a) a ~10% calibra- 
tion difference for weak sources between NVSS 
and ATCA flux densities, combined with (b) a 
small number of sources at 5 < —40° which have 
anomalously large values of Qia, perhaps simply 
because of small number statistics. 

It remains unclear whether the likely '^10% cal- 
ibration difference between Northern and South- 
ern hemispheres noted above for weak sources is 
a property of NVSS specifically or all VLA obser- 
vations. Furthermore the 4.8 GHz measurements 
of galaxies made with th e Bonn 100m single dis h 



give 94.8 = 2.71 ± 0.03 (jWunderlich et all 119871 ). 



assuming a spectral index of 0.75 this then gives 
an estimated qiA — 2.31 ± 0.03 close to the VLA 
value, supporting the idea that Northern hemi- 
sphere qiA values are consistent with each other. 
Given the importance of sorting out the origin of 
radio flux scale differences between hemispheres 
it seems that dedicated simultaneous ATCA/ VLA 
observations of sources with different flux densities 
should be made. 

4.2. VLA 4.8 GHz Maps 

The results of our VLA observations are sum- 
marized in Table 3. We successfully detected and 
mapped 95 out of the 110 observed sources. The 
15 non-detections are probably large galaxies with 
a surface brightness too low to be detected with 
the VLA conflguration used. 

The images for the VLA detections are shown 
in Figure [2] and Figures 13.1-3.71 with the former 
showing sources which are also VLBI detections 
and the latter showing those that are VLBI non- 
detections. The images sho w a wide var i ety o f 
morphologies. Following Ne ff fc Hutchingsl (jl992l ) 
we have classified the sources into four different 
morphological classes: Core (C), Core and Exten- 
sions (CE), Extended with no core (E) and Unde- 
tected (U) each class containing 28, 62, 5 and 15 
sources respectively. There are two sources with 
clear ring structures, namely IRAS0509H~0508 
(NGC1819) and 23488+1949 (NGC7771), there 
are an additional four showing elongated linear 
structures typical of edge-on disks. 
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Fig. 1. — FIR to Radio correlation plots for the COLA North sample at 1.4 GHz and at 4.8 GHz. In both 
panels sources detected by VLBI are shown circled. Left: 1.4 GHz radio spectral luminosity calculated using 
NVSS tiuxes versus FIR spectral luminosity Lfir (see §3.3|) . Most of the COLA North sources are tightly 
clustered around the median qiA — 2.34 estimated for a sample of 1809 galaxies bv lYun et al. ( 2001 ) which 
is shown by a solid line. Note on average that for the sample Lfir = 1 x 10^^ WHz"-'^ corresponds to 
LiR — Right: 4.8 GHz spectral luminosity calculated using the recovered 4.8 GHz fluxes from 

our VLA maps (see i )3.ip versus Lfir- Arrows indicate 5a upper limits for the 15 sources not detected in 
th ese observations. Crosses indicate the 20 sources having no VLBI data. The median 94 § = 2.83 measured 
bv I Wunderlich et al. (1987) from a sample of 99 normal spirals is shown by a solid line. In both panels, the 
gray area defines the region of ±0.3 dex around the median and a dashed line is drawn at — Idex. 
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2. VLA 4.8 GHz naturally weighted maps for 19 of the 20 VLBI detected sources (the well known core-dominated source 

031643+4119 (3C84) is omitted). VLA images of the sources undetected by VLBI are given in Figure [3. 1-3.71 Shown are the highest 
resolution images of each source which in all cases come from A-array observations except for 15327+2340 {Arp220) and 23488+1949 
which come from BnA observations. In 05336+5407 two features are detected 30' apart, only the brighter more compact NW feature 
which contains the VLBI source is shown. In all of the figures the IRAS source name is indicated in the top left-corner. Peak brightness 
in mjybeam"^ is indicated in the top right corner. Black contours are drawn at 10, 20, 30 and 40% of the peak brightness. White 
contours are drawn at 50, 60, 70, 80 and 90% of the peak brightness. Negative contours are drawn with dashed lines at —10, —20, —30 
and —40% of the peak brightness. A cross indicates the position of the VLBI core detections for those source having observations using 
the Eb-Ar baseline. Axes tick-marks arc in arcscconds. The horizontal line on the bottom-right corner indicates 500 pc at the distance 
of the source. 



Total recovered fluxes at 4.8 GHz (54.8 total) 
were obtained by integrating the flux over a region 
enclosed by the lowest reliable contour of radio 



emission. The median g'4.g — log(LFiR/^4.8 tot) 
obtained for the sample is 3.03 ±0.05 which is sig- 
nificantly higher than the 2.71 ±0.03 measured by 
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Fig. 3.1-3.7. — VLA 4.8 GHz naturally weighted maps for the VLA detected sources of the COLA North 
sample without VLBI detections. Maps for the VLBI detections are given in Figure [21 Column 2 in Table 3 
gives the VLA array used, for sources observed in two arrays the higher resolution image is shown. IRAS 
names are indicated on the top- left corner of every panel. Peak brightness in mJy beam~^ is indicated on the 
top-right corner. Tick marks are separated by 1 arcsecond with the restoring beam shown at the bottom-left 
corner. The horizontal line at the bottom-right corner is 500 pc at the distance of the source. Only the first 
six images are shown here the remaining images can be found in the online version Figures 3.2 to 3.7. 



Wunderlich et al.l (|l987l ) using single dish obser- 
vations of a sample of 99 normal spirals. We con- 
clude that we are on average resolving out ~50% 
of the 4.8 GHz flux density. We note however that 
for those sources with compact central features the 
fraction of missing flux density is < 20%. 

4.3. VLBI Detections 

We detect, on at least one VLBI baseline, 20 out 
of 90 sources observed. The measured flux densi- 
ties on each baseline are given in Table ID On the 
shortest Eb-Wb baseline (see Figure the signal 
to noise ratio for most detections is generally sig- 
nificantly larger than the 6ct detection threshold 
(see Figure |4]) while none detections are bounded 
to <4cr. It is notable that we obtain VLBI detec- 
tions for 5 of the 6 sources with IR luminosities 
above lO^-'^ ^'i© . A histogram showing the VLBI 
detections and non-detections versus IR luminos- 
ity is presented in Figure [5]-top which shows that 
the probability of detecting a VLBI core increases 
with IR luminosity. 



Table |4] shows that there are 9 sources with 
simultaneous measurements on the Eb-Wb and 
Eb-Jbl baselines. Of these, 7 have similar fluxes 
(within the noise) on these two baselines, the lat- 
ter of which is ~2.6 times longer than the first, 
showing that the VLBI component is unresolved. 

Of the 14 sources with simultaneous Eb-Wb 
and Eb-Ar observations 9 are detected on the Eb- 
Ar baseline, but with significantly reduced flux 
density, indicating that the VLBI sources are re- 
solved on the longest baseline. We estimate an 
angular size (see Table |4]) by approximating the 
flux ratio between both baseli nes to that exp ected 
from a circular Gaussian (see PearsonI 19991 ). 



4.4. FIR/radio ratios of the VLBI Detec- 
tions 

Figured} left shows that the total 1.4 GHz VLA 
flux densities of our VLBI detected sources fol- 
low the FIR to radio correlation. There are 
15/20 detections within ±0.3 dex o f the empiri- 
cal line defined by lYun et al.l (|200ll) ( i.e. within 
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Table 3 

Radio and optical observations and results 



IRAS VLA VLA S4.S total 94.8 S4.8 oomp ^Mxe^ log Tfc VLBI Optical 

Name Array Morph (mjy) (mJy) ["x"] [K] Epoch Class 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



00005+2140 


BnA-A 


CE 


5.66 


2.97 


0.00 






1* 


S 


00073+2538 


BnA-A 


C 


18.40 


2.85 


1.23 


0.45x0.22 


2.63 


1 


C 


00506+7248 


A 


CE 


31.60 


2.95 


13.60 


1.27x0.66 


2.75 


2 




00521+2858 


BnA 


E 


6.07 


3.13 


3.40 


3.15x0.91 


1.61 


1 




00548+4331 


A 


CE 


13.21 


2.95 


13.11 


1.51x0.41 


2.87 


2 




00555+7614 


A 


U 












2 




01503+1227 


BnA 


CE 


6.23 


3.18 


5.25 


1.58x1.03 


2.05 


1 


c 


01519+3640 


A 


CE 


4.38 


3.23 


3.92 


1.13x0.74 


2.21 


1-2 


c 


01555+0250 


BnA 


CE 


8.23 


2.92 


1.85 


2.04x1.06 


1.47 




L 


01556+2507 


BnA 


CE 


11.16 


2.95 


7.37 


1.32x1.02 


2.28 


1 


H 


01579+5015 


A 


CE 


2.49 


3.39 


1.82 


0.59x0.38 


2.45 


2 


L 


02071+3857 


A 


E 


12.20 


3.15 


0.93 


0.85x0.56 


1.83 


2 




02080+3725 


A 


U 












1-2 


H 


02152+1418 


BnA 


CE 


2.39 


3.82 


1.34 


1.49x0.94 


1.52 


1-2 




02208+4744 


A 


CE 


10.92 


3.08 


10.61 


1.78x1.02 


2.31 


2 




02253+1922 


BnA 


U 












1-2 


H 



Columns — (1): IRAS Name (2): VLA array(s) used (3): Morphology class of VLA 4.8GHz image. C-Compact, E- 
Extcndcd, CE-Compact plus Extended, CER-Compact plus extended with ring, U-Undetected. (4): Recovered 4.8GHz 
flux density in VLA image. (5). Log of ratio of 4.SGHz recovered VLA flux density to FIR flux density. (6): Total 
flux density of gaussian fitted to brightest feature in VLA image (after removing VLBI component, hence zero indicates 
that brightest feature was dominated by the VLBI component). (7): Fitted gaussian major and minor axes to peak 
VLA feature (FWHM arcsec). Note for IRAS15327+2340 (Arp 220) these are for the western nucleus as taken from the 
literature (see Section l4.5t . (7): Fitted gaussian major and minor axes to peak VLA feature (FWHM arcsec). (8): Peak 
4.8GHz brightness temperature of gaussian component (see Section 14. 5t . Note if formally unresolved along one axis a 
50mas FWHM was assumed along that axis). (9): VLBI Epoch(s) observed (asterisk indicates VLBI detection, see Tablel4] 
for details). (10): Optical Class (see Section [176l l with H=HII, S=Soyfort, L=LINER and C=Composite. 



the grey region shown). Of the remainder two 
are just shghtly above and two shghtly below the 
±0.3 dex hne. Only one source, 3C84, is signifi- 
cantly above the correlation. At 4.8 GHz a sim- 
ilar trend is observed with 14/20 detections lo- 
cat ed within ±0.3 dex of the empirical line defined 
by IWunderlich et "ah ( 19871 ). Consistent with the 
above results we find the median gi.4 for the VLBI 
detections and non-detections to be indistinguish- 
able (i.e. 2.30 ± 0.08 and 2.31 ± 0.02 respectively). 
Likewise if we omit 3C84 the mean of the VLBI 
detection is 2.33 ± 0.08 and of the non-detections 
is 2.32±0.02. We can conclude that the data show 
that excepting 3C84 the total radio flux density is 
dominated by star-formation regardless of whether 
or not a VLBI core is present. 

The above results appear to differ from those in 
COLA-S where a difference in median (71.4 between 
VLBI detections and non-detections was claimed. 



are calculated from the internal dispersion about 
th e median qi .4 for the two sub-samples as given 



From Table-3 of lCorbett et all (|2002l ) the median 
(7i.4 value and its estimated one sigma error for 
the nine definite COLA-S VLBI detections was 
2.07 ± 0.17 while for the VLBI non-detections it 
was 2.44 ± 0.03. The errors on theses estimates 



Corbett et al. 



(|2002l ). We ascribe this appar- 
ent two sigma difference between the two sub- 
samples in the South as being due to the lower 
VLBI sensitivity in the South combined with the 
overall form of the gi.4 distribution. This distri- 
bution has a tail of radio core-dominated, low gi.4 
sources (such as 3C84 in the North and NGC5793 
in the South) which are easily VLBI detected plus 
weaker VLBI cores in galaxies whose radio emis- 
sion is star-formation dominated. The lower sen- 
sitivity VLBI observations in the South contain a 
larger fraction of the former class of source which 
reduces the estimated (71.4 of VLBI detection com- 
pared to non-detections. Specifically in the South 
4 out of 9 VLBI detected sources have (71.4 < 2 
while the fraction is only 3 out of 20 in the North. 
Because low gi.4 sources comprise almost half of 
the VLBI detected sample in the South they can 
bias the median (71.4 while in contrast in the North 
such outlier sources have negligible effect. 
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Table 4 
VLBI Results 



COLA 


Other 


Eb-Wb 


Eb-Jbl 


Eb-Ar 


Angular Size 


Linear Size 


Optieal 


name 


name 


[mJy] 


[mJy] 


[mJy] 


[mas] 


[pc] 


Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


00005+2140 


Mrk 334 


1.67 (0.15) 




1.02 (0.03) 


0.68 


0.29 


s 


02568+3637 




1.61 (0.14) 


1.92 (0.16) 








s 


03164+4119 


3C84 






15500(0.02) 






L 


04435+1822 




3.32 (0.14) 




3.10 (0.03) 


0.25 


0.07 


H 


05054+1718 




1.64 (0.18) 




< 0.20 


>1.41 


>0.50 


S 


06140+8220 




4.12 (0.32) 












06239+7428 






1.51 (0.19) 










F07258+3357 


N2388 


1.38 (0.18) 


1.26 (0.14) 


< 0.20 


>1.35 


>0.35 


H 


08354+2555 


Arp 243 


2.03 (0.10) 




0.27 (0.02) 


1.38 


0.49 




13126+2452 




1.38 (0.17) 




0.88 (0.04) 


0.65 


0.16 




13183+3423 


Arp 193 


2.17 (0.17) 


3.85 (0.15) 


0.70 (0.04) 


1.03 


0.46 


C 


13470+3530 




1.79 (0.11) 


5.64 (0.16) 


< 0.18 


>1.47 


>0.48 




15327+2340 


Arp 220 


4.37 (0.17) 




2.68 (0.04) 
1.02 (0.04) 


0.40 


0.14 


L 


16478+6303 


N6247 


2.60 (0.12) 


< 0.82 










16577+5900 


Arp 293 


6.11 (0.27) 


5.36 (0.17) 








C 


17468+1320 




8.33 (0.21) 


5.31 (0.13) 


< 0.25 


>1.82 


>0.57 




17548+2401 




4.45 (0.20) 




0.39 (0.04) 


1.52 


0.58 


L 


23007+0836 


N7469 


5.12 (0.22) 


6.32 (0.15) 








S 


23488+1949 


N7771 


< 1.14 




0.33 (0.04) 


<1.08 


<0.28 


L 


23488+2018 


Mrk 331 


5.02 (0.21) 


4.92 (0.13) 


0.30 (0.04) 


1.63 


0.57 


C 



Note. — Columns — (1): IRAS Name (from Point Source Calalog, excepting F0725S+3357 which is taken from 
the Faint Source Catlog). (2): Other name. (3)-(5): VLBI flux and la uncertainty (in brackets) obtained from the 
peak in the Delay— Rate map. Upper limits for the non-detections are indicated at 5cr. A dash is shown for sources 
not observed at that baseline. For IRAS15327-f 2340 (Arp 220) on the longer baselines the two rows give results 
for the West and East nuclei respectively. (6): FWHM of a circular Gaussian source fitting the fiux ratio between 
the Eb— Wb and Eb— Ar baselines (see ^4.3| l . Sources with flux upper limits on the Eb— Ar have consequently size 
lower limits. (7): Corresponding fltted Gaussian FWHM linear sizes. (8): Optical Class (see Section |4.6| l with 
H— HIT, S— Scyfert, L— LINER and C— Composite; a dash indicates no optical line data. 
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Fig. 4. — Radio spectral luminosity at 4.8 GHz on the Eb-Wb baseline for VLBI detected sources (excluding 
3C84) versus IR luminosity. The source Arp 220 is shown as a filled circle. The scale on the upper horizontal 
axis indicates the core-cq l lapse supernova rate vsn estimated from the IR luminosity using the formula in 
Ivan Buren fc Greenhous j ( 1994[) (after applying a small FIR to IR luminosity correction). The luminosity 
at maximum light measured for the powerful Type Iln SN1986J is indicated with a horizontal dashed line at 
2.1x10^^ WHz~^. The diagonal line represents the spectral luminosity detected on the Eb-Wb baseline for 
Arp 220 scaled by IR luminosity. The lower horizontal dotted line indicates our typical 6(T=780 fxJy detection 
threshold on this baseline at the median distance of the COLA sample. 



4.5. Relation Between VLA and VLBI 
Sources 

Visual inspection of the VLA radio maps indi- 
cates a bias for VLBI detections to be in sources 
with compact structure and high peak brightness 
(see Figure [2]). This is confirmed in terms of VLBI 
detection rate versus morphology class (see §4.2p 
. Amongst the VLBI observed sources there are 
13/27 detections of class C, 6/48 of class GE and 
1/5 of class E. We do not detect any class U source. 
This result was confirmed by plotting VLA visibil- 
ity amplitude versus uv distance. The majority of 
VLBI detections were indeed found in sources with 
compact, but not unresolved, arcsec-scale struc- 
tures. The size of these structures was estimated 
by subtracting an unresolved source equivalent to 
the VLBI detection from the brightest component 
in each VLA map and then fitting a Gaussian to 
the difference. 



In two cases (00005-^2140 and 06239-^7428) the 
VLBI core subtraction removed all the flux den- 
sity implying that the brightest VLA component 
was caused completely by the VLBI source it- 
self. We consider these two sources together with 
03164-1-4119 (3G84) to be VLBI core- dominated 
sources and omit them from the subsequent anal- 
ysis. For the rest of the sources the task JMFIT 
was used to fit a Gaussian within a tight box con- 
taining the remaining feature. For the vast ma- 
jority of sources in classes G and GE with a clear 
compact component there was no ambiguity in the 
Gaussian fitting, i.e., the results of the fits were the 
same regardless of the size of the fitting box. For 
sources in class E the size estimates were found 
to vary by up to a factor of ^2 depending on the 
chosen box. However in all such cases the fitted 
sizes were large (^2"). The resulting component 
flux densities and major and minor axes 0m and 
9m are listed in columns 6 and 7 of Table 3. 
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log([NII]/Ha ) log([SII]/Ha ) log([OI]/Ha ) 



Fig. 6. — Optical-line diagnostic diagrams ("BPT" diagrams) for the COLA-N som'ces that are both VLBl 
observed and have available optical line data. Dots indicate VLA detected sources and crosses those un- 
detected by the VLA. VLBI detected sources are show n circled, a) [N11]/H q! v/s [OlllJ/H/3 diagram. The 
solid and dashed c urves are respectively the theoretical iKewlev et al.l ( 20011 ) extreme starburst line and the 
iKauffmann e t aD ( 2003 ) empirical line for star- forming galaxies. Sources that lie between the two curves 
are classified as composite Starburst/AGN. b) [SII]/Ha v/s [0III]/H;3 diagram c) [OI]/Ha y/s \Olll]/ll0 
diagra m. In the latter two panels the solid curved line is the extrem e starburst line fro m iKewlev et al 



( 2001 ) and the straight diagonal line separates Seyferts from LINERS ( Kewlev et al. I [20061 ). Note that for 
some sources not all line ratios are plotted because some lines were not detected. The shaded VLBI detected 
source is 15327-f 2340 (Arp220). 



Peak brightness temperatures Th at 4.8 GHz for 
the compact VLA structures were calculated using 
the relation Tb ~ 35x(54.8peak)(6'Mx6'™)-i K° 
(Columns in Table 3). In a few sources the fitted 
deconvolved minor axis was zero, in these cases 
its size was set to 50 mas which was equal to the 
smallest size found in other sources. 

Uniquely amongst sources with compact struc- 
ture we did not collect any VLA A-array data 
for the well known ULIRG source 15327-^2340 
(Arp220). In order to treat this source in a simi- 
lar way to the rest we give in Table 3 radio compo- 
nent properties from the lite rature for the b r ighter 
western nucleus (sizes from ICondon et al.l (|l991l ) 
with flux de nsities scaled t o 5 GH z using the spec- 



tral index of Rovilos et al. Ijiooi). 



As can be seen from Figure [5] (bottom panel) 
VLBI detections seem to be strongly correlated 
with those sources whose VLA fitted components 
(after removing any VLBI contribution) have the 
highest Tf). VLBI detectability is less strongly cor- 
related with large Lir (top panel). Two sided 
Kolmogorov-Smirnov tests comparing the distri- 



bution of VLBI detected and non-detected sources 
in LiR and T;, in each case rejects the hypothesis 
that VLBI/non-VLBI sources are drawn from the 
same population at probability levels of P = 0.82 
and P — 0.98 respectively. This result demon- 
strates that the stronger correlation for VLBI de- 
tectability is with Tb rather than with Lib.. 

4.6. Optical Spectroscopy Classifications 

To classify sources based on their line ratios 
from our optical spectrosc opy observations (se e 
gSH) we used the method of IKewlev et~all (|2006f ). 
This scheme first classifies sources into HII galax- 
ies, AGN/star- formation composites and pure 
AGN bas ed on a mixture of t heoretical (maximal 
starburst, IKewlev et al.l 2001 ) and empi rical line- 
ratio boundaries ( Kauffmann et al ]|2003l) delineat- 
ing pure star-forming galaxies. These boundaries 
are shown respectively as solid and dashed curves 
delimiting the bottom left corners of the diagnos- 
tic line ratio diagrams in Figure|6l HII galaxies 
(class H) are to the left and below all the curves 
including the dashed line in Figure [S^. Composite 
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Fig. 5. — Histograms of IR luminosity (top panel) 
and peak brightness temperature of VLA scale 
radio component (bottom panel) for VLBI ob- 
served sources. In both panels VLBI detections 
are shown black. VLBI non-detections for sources 
respectively detected and undetected by the VLA 
are shown white and grey respectively. The bright- 
ness temperature in the bottom panel is calculated 
after subtraction of the VLBI component from the 
VLA map (see ! i4.2|) . Note for ease of compari- 
son the three VLBI core-dominated detections (i.e 
00005-^2140, 03164-^4119, 06239-^7428) are not 
plotted in either panel. A much stronger correla- 
tion between VLBI detection and high brightness 
temperature VLA scale emission is found than 
with IR luminosity (see Section 14.51 and also Fig |8] 



sources (class C) lie between the solid and dashed 
lines in the [Nil] /Ha diagram and AGN are above 
and to the right of the curves in all diagrams. 
The AGN population is further divided between 
Seyfert (class S) or LINER (class L) using the 
[SII]/Ha and [01] /Ha diagnostic di a grams and 
the empirical results of iKewlev et all f|2006h who 
found two well defined populations tracing differ- 
ent position angles on these diagrams. For a signif- 
icant number of our sources either the [SII] or [01] 
lines were not detected. In order to increase the 
number of detections we therefore us ed the same 
2- of- 3 criterion for classification as lYuan et al 



( 2010l ) whereby if there was no data for a given 
line ratio or if one diagram gave a classification 
inconsistent with the other two we chose the clas- 
sification given by the majority of the diagrams. 

Our available optical data allowed for the classi- 
fication of 66 sources as summarized in column 10 
of Table [31 The fraction of different classifications 
agrees well with that found amongst s ources of 
similar luminosity by f Yuan et al. ( 2010t ) . For the 
subset of sources observed by VLBI the compact 
radio core detection rates for each optical cate- 
gory are; Seyferts 4/8 (50% ), LINERS 4/6 (66 %), 
Composite 3/20 (15%) and HH 2/15 (14%). These 
VLBI detection fractions are broadly consistent 
with the idea that VLBI cores trace AGN activ- 
ity since the detection fraction is largest amongst 
AGN optical types (see 

5. Discussion 

5.1. Nature of Detected VLBI Sources 

The primary objective of our VLBI observa- 
tions was to detect compact cores indicating AGN 
activity. Here we discuss whether our 20 VLBI 
detections can be taken as unambiguous signs of 
AGN activity or whether they might be confused 
by the high spatial frequency tail of arcsecond 
scale diffuse emission, radio supernovae (SNe), or 
supernova remnants (SNRs). 

5.1.1. VLBI Brightness Temperature and Diffuse 
Starburst Emission 

Could the VLBI detections be due to spatial 
substructure within the general ISM starburst- 
powered radio emission? As noted in §4.51 the fiux 
densities on the Eb-Wb and Eb-Jbl baselines are 
usually comparable so implying a source size less 
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than about half of the Eb-Jbl fringe spacing, i.e. 
<10mas. A typical VLBI detection flux density 
of ~3 mjy would then imply a brightness tem- 
perature of at least 10^'^ K. Even larger bright- 
ness temperatures are found on the Eb-Ar base- 
line where sub-milliarcsecond sizes imply bright- 
ness temperatures of >10®K. These high bright- 
ness temperatures are inconsistent with star- 
formation powered synchrotron emission which 
is predi cted to have bright ness temperatures of 
<10^K (|Condon et al.lll99ll) . 

Given that we only have fringe detections on 
single baselines it is possible that such detections 
(especially on the Eb-Jbl baseline) could be due to 
star-formation powered radio emission combined 
with an unusual source brightness distribution. 
One such distribution would be a thin rectangle 
with its small dimension < 10 mas along the di- 
rection of the fringe spacing; this would however 
require both unusual brightness distributions and 
chance alignments in position angle in a number of 
sources to affect our results. A second such distri- 
bution would be one with order 100% amplitude 
fluctuations of lO'^ K brightness on scales < 10 mas 
over a region of order 1000 mas in diameter, giving 
large r.m.s visibilities of the same size on Eb-Wb 
and Eb-Jbl. The ratio of detected visibilities on 
these two baselines seems too close to unity to be 
explained by such a random model. Additionally 
such a brightness distribution seems very unlikely, 
one would instead expect larger fluctuations as we 
go to larger separations within the source giving 
rise to larger r.m.s visibilities on the shorter base- 
line. This is what is in fact observed in Arp 220. 
Taking the data of iParra et al.l (|2007l ) and remov- 
ing the effect of the compact sources we detect the 
diffuse emission on Eb-Wb but not on the Eb-Jbl 
baseline. Since Arp 220 is brighter and more com- 
pact than other COLA sources it seems unlikely 
that the Eb-Jbl detections in these can be due to 
randomly fluctuating diffuse lO'^ K emission. 

5.1.2. Compact Sources in Arp 220 

For the case of Arp 220 (15327+2340) most 
of the flux density detected in our VLBI snap- 
shot observations clearly comes from non-AGN 
activit y. This galaxy has been known for many 
vears (ISmith et al.l Il998l : iLonsdale etUI 120061 



of luminous radio supernova remnants (SNRs, in- 
teracting with the ISM) and supernovae (SNe, 
interacting with progenitor winds). Arp 220 is 
however a very unusual source within the COLA 
sample being by far its most IR luminous object. 
It follows that a confirmed SNR/SNe origin for the 
VLBI snapshot fringes in Arp 220 does not nec- 
essarily imply that the same explanation is valid 
for the rest of the VLBI detections. The diag- 
onal line in Figure |3] shows the Eb-Wb spectral 
luminosity of Arp 220 scaled linearly with Ljr and 
demonstrates that if the VLBI contribution from 
SNR/SNe scales linearly with star-formation rate 
then these cannot explain the VLBI detections. 
In the next subsections we discuss in more detail 
whether SNRs and SNe can in fact explain the 
VLBI detections. 

5.1.3. Supernova Remnant Origin 

Chomiuk fc WilcotsI (|2009t ) find in their study 



Parra et al.l 2007 ) to harbor multiple compact ra- 
dio sources consistent with a mixed population 



of the radio luminosity function of SNRs in nearby 
normal and starburst galaxies a constant power 
law exponent and a scaling factor proportional 
to SFR. They also find both theoretically and 
empirically that the luminosity of the brightest 
SNR in each galaxy scales linearly with the SFR. 
Since even on our shortest VLBI baseline bright 
SNRs will be separated by several fringe spacings 
their expected contribution to the correlated VLBI 
flux density will be an incoherent sum of all the 
SNR flux densities. This incoherent sum will how- 
ever, because of the c onstant slope SNR luminos- 
ity function found bv IChomiuk k. WilcotsI ( 2009h 
simply be some multiple of the brightest SNR flux 
density. Since this latter quantity is linearly corre- 
lated with SFR we expect the total observed VLBI 
signal due to SNRs to also linearly scale with SFR. 
Since about half of the flux density seen on the Eb- 
Wb baseline for Arp 220 comes from SNR a VLBI 
contribution which lies a factor of two below the 
diagonal line shown in Figure |4| is therefore ex- 
pected; in all cases this lies below the luminosity 
of the detected VLBI cores. 

In fact estimating the SNR contribution to 
VLBI flux density in COLA galaxies by scal- 
ing Arp 220 by Lm as described above proba- 
bly results in an overestimate of the SNR con- 
tribution for COLA galaxies. This is because 
IChomiuk fc WilcotsI (j2009) find that ahhough the 
SNR luminosity function in Arp 220 has the same 
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slope as in other galaxies the scaling factor be- 
tween SNR brightness and SFR is larger (pos- 
sibly because of the very high gas density in 
Arp220). The true VLBI SNR contribution in 
COLA galaxies although having the same linear 
slope as drawn in Figure|4] therefore likely passes 
below the Arp220 point; because this source is 
over-luminous in SNR for its SFR. 

5.1.4- Radio Supernova Origin 

A difficulty for calculating the possible contri- 
bution of radio SNe to our detected VLBI flux 
density is that unlike the case for radio SNR no 
well defined radio SNe luminosity functions have 
been published. It is clear however that only 
the mo st luminous Ty peic or Typelln optical 
classes (jChevalieii 120061 ) would be detectable at 
the distance of COLA galaxies. The most power- 
ful Typeic objects are often associated with 7-ray 
bursts and have rapid rise and decay times. The 
slightly weaker Type Iln SNe are much longer lived 
and hence more likely to be detected. Optically 
about 4% of core-c ollapse SNe are of Typelln 
( Smartt et al. 2009[ ) however even these show a 
wide range in radio luminosity. For instance 
all of the six Iln's observed by Ivan Dvk et al 



( 19961 ) were undetected in the radio with 4.8 GHz 
spectral luminosity <1.4xl0 ^^ WHz~^ The 



likely Type Iln's in Arp220 (jParra et all 120071 ) 
have spectral luminosity of ~ 4x10^° WHz"^. 
The most radio luminous Typelln so far de- 
tected (2x1 0^1 WHz-i at 4.8 GHz) is SN1986J 



in NGC891 (jRupen et al.lll987[ ). Sources close to 
this maximum seem however to be rare. Consid- 
ering the number of opti cally identified Type Iln's 
searched in the radio (see Chandra et al.lliooii and 
references therein) in normal galactic disks we can 
estimate that at most one in ten Type Iln's (or 
0.4% of all core collapse SNe) reaches to within a 
factor of two of the SN1986J radio luminosity. 

Amongst COLA VLBI detections there are 
six objects other than Arp220 with luminosities 
above that of SN1986J at maximum light (i.e. 
above the dashed line in Figure|4]). A single RSNe 
origin seems to be ruled out for these objects and 
we are forced to consider models consisting of 
two or more extremely luminous objects. Such 
SN1986J-like objects stay within a factor of two 
of their peak 4.8 GHz luminosity for ~5 years and 
host galaxies of these VLBI detections have to- 



tal supernova rates of order 0.5year~-^. If 0.4% 
of core collapse are SN1986J-like, then we expect 
0.01 such objects on average per COLA galaxy. 
The probability of getting just one example of two 
or more such luminous objects in a galaxy within 
our whole sample is thus only 1.5%. 

FigureHshows 7 VLBI detections at 1 x 10^° W Hz~i 
that could potentially be explained as single 
SN1986J-like objects observed when above 50% 
of maximum light. Assuming the same radio lu- 
minous SNe lifetimes and fractions as before and 
typical SFRs of this group of ^0.25year~^ there 
should only be 0.5 such objects in the whole sam- 
ple. Finally there are three VLBI detections at 
~ 0.6 x 10^° W Hz~^ , it is possible that one of these 
objects could be a SNe, though without a detailed 
radio SN luminosity function it is hard to make 
any definite statements. 

A general argument that suggests that the vast 
majority of our VLBI detections are not explained 
by SN comes from the fact (see that all 

our VLBI detections are well above our detection 
limit. If the origin of the VLBI detections were 
single or multiple RSNe such a gap in luminosity 
between our weakest detection and the sensitivity 
limit is hard to explain. A single RSNe spends 
a much longer period at low luminosities than at 
high, and intrinsically weaker SNe are more com- 
mon than strong ones, so one would expect many 
more detections just above the detection cutoff 
than are observed. 

We conclude that, while one or two of the less 
luminous detections could be caused by radio SNe, 
the vast majority of our VLBI detections cannot 
and must instead be AGN powered. 

5.2. Nature of Extended Radio emission 

For most VLBI detected COLA sources (17/20), 
a large fraction of the total radio flux density arises 
from several 100 pc scale structures seen in the 
VLA maps. As noted in i)4.1l the total flux den- 
sities of nearly all COLA sources, including both 
VLBI detections and non-detections (see tj4.4l and 
Figure[T|) , closely follows the FIR to radio correla- 
tion at both 1.4 and 4.8 GHz. This result indicates 
that most of the radio emission recovered in the 
VLA maps is powered by star formation activity. 
It is known from observations of the starbursts 
in M82 and Arp 220 that the associated star for- 
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mation induced compact emission from SNR a nd 
SNe is just a few percent (ILonsdale et alJl2006l ) of 
the total diffuse component radio emission; hence 
there is no contradiction between star-formation 
being the main power source for the total radio 
emission and our conclusion that it does not pro- 
vide enough correlated flux density to explain the 
VLBI detections. 

5.3. Nature of integrated IR emission 

Figure[7] shows a IR color-color diagram for the 
90 sources with VLBI observations in the COLA 
North sample. Any position on this diagram can 
be reached by a linear combinatio n of starburst 
temp erature, AGN and reddening ( Kewlev et al.l 
2000). It is therefore impossible to unambiguously 
estimate the proportions in which these ingredi- 
ents are combined. However, most of the sources 
in the sample appear to follow the line traced by 
blackbodies of progressively increasing tempera- 
tures (i.e. the starburst line) with a scatter ap- 
parently increasing with temperature suggesting a 
larger AGN contribution to the IR energy bud- 
get in warm starbursts. The overall scatter of 
the GOLA samp le is howeve r much smaller than 
that observed bv lRush et al] (|l993l ) in their sam- 
ple of optically selected Seyferts. We therefore ar- 
gue that for most sources in our sample that their 
IR emission is mostly powered by starburst activ- 
ity with any AGN contribution being energetically 
unimportant. 

All 10 VLA non-detections (see g^J with VLBI 
data are clustered near the cold end of the star- 
burst line (crosses in Figure[7]). As discussed in 

2 1 these are likely to be large sources whose rel- 
atively large FIR luminosity is not due to con- 
centrated nuclear starburst activity but is instead 
caused by the aggregation of normal star forming 
regions spread over their large galactic disks. 

5.4. Comparison of optical and radio source 
classifications 



As discussed in ij4.6l we find that the fraction 
of VLBI detections in AGN-like optical types is 
substantially higher than in HII type galaxies, 
supporting the interpretation that VLBI detec- 
tions indicate AGN activity. We find however 
no one-to-one correspondence between the opti- 
cal classification and VLBI detection. The fact 



that half of the observed Seyfert sources (4/8) 
w ere not detect e d by VLBI supports the result 
of lGorbett et al.l (|2003l ) who found observationally 
two classes of Seyferts with and without detectable 
VLBI radio cores. This could simply be an artifact 
of our limited VLBI sensitivity or it could mean 
that there are physically two classes of Seyfert 
with different degrees of radio-loudness. The rela- 
tively high VLBI detection rate of LINER sources 
(4/6) is consistent with m odels in which this c lass 
is usually AGN powered ( Kewlev et al. 200 6^ al- 
though it is clear that in some sources LINER- 
like spectra can arise in shocks , i.e. for the case 
of Arp220 (|Arribas et al.ll200l[ ). The relatively 



low incidence of VLBI detections in HII galaxies 
(2/15) suggests that cases of completely optically 
obscured AGN are rare. Combining the VLBI 
and optical results we find the overall incidence of 
AGN activity to be high. If we add the four non- 
VLBI detected Seyferts to the 19 VLBI detections 
(excluding Arp 220) we get a minimum fraction of 
sources containing AGN of 22/90 (25%). 

5.5. Implications for the AGN starburst 
connection 



An interesting correlation found in the present 
work is that between the presence of a VLBI de- 
tection and of high brightness arcsecond-scale ra- 
dio emission (see ij4.5l and Figure[S]) . In H5.1\ we 
argued that VLBI detections are due to AGN ac- 
tivity rather than SNe or SNRs. In i )5.2l we further 
argued that the extended radio emission seen with 
the VLA is due to star formation. Taken together 
these results imply a strong connection between 
AGN activity and the presence of a compact star- 
burst. 

A significant body of evidence in favor of a 
strong starburst / AGN connection has been accu- 
mulating in recent years. In nearby AGN evidence 
is seen for for nuclear post-st arburst stellar pop- 
ulations of age ^10*^-10^ year (|Davies et al.ll2007l ) 
suggesting a possible evolutionary connection with 
nuc lear starbursts preceding AGN activity (see 
also Heckmanll2008[ ) . There are however also many 
individual e xamples of concur rent starburst/ AGN 
activity (e.g. Fathi et al1l2006 ). The present paper 
has found evidence linking AGN activity with con- 
current compact starbursts in FIR selected sources 
and it is interesting to consider possible physical 
interpretations of this result. 
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Fig. 8. — Starburst FIR equivalent brightness Sfir(4.8) versus central starburst component radius R for 
sources both imaged by the VLA and observed by VLBl. VLBI detections are in dicated by a circ l e with 



Arp220 emphasized as a filled circle. This diagram duplicates a similar one in [Thompson et al.l ( 20051 ) 
for more luminous IR galaxies. In this plot Sfir(4 g) = iFiR(4.8)/^radio where iFiR(4 g) is the equivalent 
starburst FIR spectral luminosity estimated from the 4.8 GHz data and ^radio is the area of the central 
starburst source. Lfir(4.8) is calculated by transforming the 4.8 GHz spectral luminosity of the brightest 



VLA component L4 



Speak 



to an equivalent FIR spectral luminosity using the FIR-radio correlation. 



The Eddington-limit ed starbursts discussed by 



Thompson et al.l (j2005l ). which can also feed cen- 
tral AGN, are predicted to have central peak 
bolometric brightness of 10^'^ Lq kpc~^ consis- 
tent with the lumin ous compact starbursts from 
Condon et al.l (|l99ll ). Figure|5] plots the estimated 
FIR surface brightness versus fitted size for the 
arcsecond-scale radio structures for our COLA 
source s, matching Figure 4 of Thompson et al.l 
POOSI ) . We find that our VLBI sources are mostly 
detected in sources with bolometric brightness 
approaching or around lO^'^ Lq kpc~^ which coin- 
cides with values predicted for Eddington-limited 
starbursts. 



Thompson et al.l ( 2005 ) also argued for two 
separate modes by which such compact star- 
bursts can feed a central black hole, supply- 
ing either 0.1M©year~^ or 4M0year~^ depend- 
ing on whether the gas supply rate at several 
hundred parsec radius is smaller or larger than 
22OM0year~^. The rest of this mass is converted 
into stars. Since the star formation rates of COLA 
sources are ^220 M^gyear"^ the lower black hole 
supply rate applies. Assuming an efficiency of 10% 
this supply rate translates into an AGN bolomet- 
ric luminosity of l.SxlO^^L© which is comparable 
to the IR luminosities of the COLA galaxies with 
VLBI detections. Such a high AGN luminosity, 
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Fig. 7. — IR color-color diagram for COLA North 
sources observed by VLBI. VLBI detections are in- 
dicated with open circles and VLA non detections 
as crosses. The gray line is traced by blackbod- 
ies of temperatures progressing from 60 to 80 K 
from the left. The black line is a reddening line 
starting from the average warm Seyfert 1 nucleus 
(indicated by a star) . Squares over t his line are lo- 
cated at T25^an=l, 2, 3, 4 and 5 (see lDopita et al 
Il998l) . 



which would mostly emerge in the IR, is inconsis- 
tent with the fact that COLA sources lie on the 
FIR-radio correlation and have IR colors typically 
found in starburst powered sources (see 15.31) . We 
can speculate that a constant fraction of gas sup- 
plied to the black hole might be lost to winds or 
jets formed close to the black hole event horizon, 
thus reducing the AGN bolometric luminosity. 
The detected core radio emiss ion would be gener- 
ated at the base of these jets ( Falcke fc BiermannI 
19951 ). The starburst regulation of the gas supply 
would ensure that these jets had similar mass out- 
flow rates and mechanical luminosities therefore 
explaining the relatively narrow range of radio 
luminosity of our VLBI detections. 

The high incidence we find of AGN activity in 
compact high brightness starbursts (i.e. VLBI 
detections in 13/18 sources with specific star- 
formation rates >10^^ Lq kpc~^, see Figure 8) sug- 
gests only a short interval between the triggering 



of such starbur sts and the onset o f AGN activ- 
ity. In contrast lOavies et al.l ( 20071 ) argue (based 
on AGN stellar population studies) that there is 
a 50-100 Myr delay between the peak of nuclear 
starburst activity and the onset of AGN activity. 
These studies were however done on AGN selected 
sources. AGNs plausibly last much longer after 
the initial feeding event than does starburst ac- 
tivity. Possibly this is because once the central 
~lpc scale accretion disk is fueled with gas via 
a compact sta rburst long viscous ac cretion times 



(up to 10" yr. iKing k Pringld 120071) allow AGN 



activity to continue long after starburst activity 
has faded. Given this timescale difference ex- 
amples of currently active nuclear starbursts in 
AGN selected samples might well be expected to 
be rare. Additionally during the initial composite 
AGN/starburst phase optical obscuration may be 
very large making such objects hard to be opti- 
cally identified as AGN. Since COLA is FIR se- 
lected and we have used radio observations as the 
primary means to detect AGN both of the above 
biases are avoided allowing the detection of the 
initial phases of composite AGN/starburst activ- 
ity. 



Bo th our own observations and those of l Condon et al 
( 199l[ ) show a strong correlation between high 
FIR luminosity and compact high bri ghtness star- 



bursts . Within the IRAS-BGS sample lCondon et al 
(fl99l found that all starbursts with Lfir > 
10^^'^ Lq were in fact of this high brightness 
Eddington-limited type. If the correlation be- 
tween AGN powered VLBI sources and high 
brightness shown in Figure[5] also holds for this 
more luminous sample then one would expect all of 
these to have AGN activity detectable via VLBI. 
It can be argued that Arp 220 is a clear counter- 
example since it has high radio brightness and 
although detected by VLBI this detection is due 
to SNe and SNRs rather tha n an AGN c ore or 
lobes. However as noted by iParra et al.l (|2007[ ) 
within the western nucleus of Arp 220 there exist 
several candidate sources for a radio AGN. Recent 
millimeter continuum and ^ ^C0(2-l) spect r al line 



inte rferometry result s bvlDownes &: EckartI (l2007|) 



and ISakamoto et al.l (|2008[ ) strongly bolster the 
case that a powerful AGN does in fact exist in the 
western nucleus of Arp 220. 

More generally within s napshot 1.4 GH z VLB I 
observations conducted by iLonsdale et al.l (|l993l ) , 
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about 50% of luminous IRAS-BGS sources were 
detected on global baselines. These detections had 
however ambiguous interpretations being explain- 
able by both AGN and radio SNe/SNR models. 
In general, at other wavebands attempts to defi- 
nitely detect or reject AGN contributions in these 
very luminous starbursts have als o often proved in- 



The 15 non-detected sources are likely ex- 
tended (resolved out) galaxies with radio surface 
brightnesses below our sensitivity (see §4.2[) . 



conclu sive. If as predicted by the [Thompson et al 



(|2005h model for Eddington-limited starbursts the 
gas supply to the central black hole is relatively 
constant for compact starburst sources we may ob- 
tain the same AGN luminosity regardless of Lir. 
It follows that in more intense higher IR luminos- 
ity starbursts it may be harder to detect the AGN 
component in the midst of all the starburst ac- 
tivity. Paradoxically it may be that within the 
moderate luminosity COLA sample a strong cor- 
relation between intense starburst activity and ra- 
dio AGNs is easier to detect than in more lumi- 
nous samples. This is firstly because COLA strad- 
dles a luminosity range were there are large num- 
bers of both low brightness and Eddington-limited 
starbursts. Secondly, Lir is not high enough (ex- 
cept for Arp 220) that there is significant confusion 
in the radio between AGN and starburst induced 
structure. 

It is clear that both the VLBI snapshot ob- 
servations of COLA-N presented in this paper 
and the o bservations of IRAS-B GS sources car- 
ried out by iLonsdale et al. ( 19931 ) suffer from lack 
of uv coverage. New long track VLBI observa- 
tion which produce good quality images would 
help distinguish star-formation origins of VLBI 
flux density (by for instance seeing multiple com- 
pact VTK[_coniponents indicating star-formation, 
e.g. IParra et al.l [2007) from AGN origins (core- 
jet morphology sources). Observations of half of 
the VLBI detected COLA-N sources have already 
been made using the EVN at 6 cm wavelength. 
These data are presently being reduced. 

6. Conclusions 

The main conclusions of this paper are: 

1 We have made 4.8 GHz VLA images for 95 of 
the 110 sources in the COLA North sample. Of 
these sources 28, 62 and 5 sources have respec- 
tively Core, Core4-Extended and Extended ra- 
dio morphologies as defined by iNeff fc Hutchingsl 



2 We argue that the properties of the 4.8 GHz 
radio emission structures imaged on scales of a 
few 100 pc are consistent with being powered by 
star formation activity because the recovered fiux 
density in these structures follows the FIR to radio 
correlation for star forming galaxies (see M.l^ . 

3 We detect 22% (20 /90) of the sources observed 
in the VLBI snapshots (see Assuming that 
these VLBI detections are single Gaussian compo- 
nents the simultaneous detection at multiple base- 
lines in most cases constrains the size of the emit- 
ting region to parsec scales. 

4 The majority of the VLBI detections are well 
above our detection threshold (see Figure |4]) sug- 
gesting a lower limit of ^10^^ WHz^^ for the lu- 
minosity of these compact sources. 

5 Based on luminosity arguments we can rule out 
a radio supernova remnant (SNR) origin for all of 
the VLBI detections (see ^5X3)) . 

6 We can also rule out a radio supernova (SN) 
origin for most VLBI detections unless we assume 
that the incidence of powerful radio supernovae is 
much larger in COLA sources than in most galac- 
tic disks. We estimate that at most one or two of 
the VLBI detections could be due to radio SN. 

7 Our VLBI detections are preferentially found 
in sources with the highest radio brightness 
arcsecond-scale radio emission, even after sub- 
tracting off the VLBI cores. Since this arcsec- 
scale radio emission is powered by star-formation 
activity (see point 2 above) we conclude that there 
is a relationship between compact starbursts and 
A GN activity. This is consistent with the proposal 
of lThompson et al. 112005) that Eddington-limited 
nuclear starbursts can feed AGN. 

8 The high AGN incidence amongst compact 
starbursts implies no significant delay between the 
period of intense nuclear starburst activity and the 
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onset of AGN activity. If the AGN activity contin- 
ues long after the starburst has faded this would 
explain the predominance of post starburst stellar 
populations in AGN sources. 
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23007+0836 


ARP298 


23 03 15.62 


08 52 26.06 


4892 


70.3 


11 


,63 


23.02 


25.32 


2.30 


23106+0603 




23 13 12.75 


06 19 18.11 


3536 


51.0 


10, 


,69 


21.52 


24.29 


2.76 


23157+0618 




23 18 16.28 


06 35 09.03 


4966 


70.9 


11 


,05 


22.50 


24.80 


2.31 


23213+0923 




23 23 54.07 


09 40 02.10 


3559 


51.1 


10, 


,58 


21.72 


24.33 


2.60 


23262+0314 


ARP216 


23 28 46.67 


03 30 41.00 


5138 


73.3 


11 


,13 


22.55 


24.83 


2.28 


23414+0014 


N7738 


23 44 02.04 


00 30 59.84 


6762 


96.3 


11 


.11 


22.57 


24.86 


2.29 


23445+2911 


ARP86 


23 47 02 


29 28 23 


4845 


69.1 


10. 


,95 


22.32 


24.71 


2.39 


23485+1952 


N7769 


23 51 03.98 


20 09 01.52 


4211 


59.9 


10. 


,99 


22.41 


24.77 


2.36 


23488+1949 


N7771 


23 51 24.85 


20 06 42.30 


4256 


58.7 


11. 


,31 


22.77 


25.08 


2.32 


23488+2018 


MRK331 


23 51 26.74 


20 35 10.33 


5541 


78.7 


11. 


,49 


22.72 


25.24 


2.52 


23591+2312 


Taffy 


00 01 41.91 


23 29 45.00 


4353 


61.8 


10, 


,91 


22.79 


24.66 


1.88 



Columns — (1): IRAS name (from Point Source Calalog, excepting F07258+3357 which is taken from the Faint Source 
Catlog). (2): Other name. (3)-(4): Right Ascencion and Declination (from VLA observations, unless undetected in which 
case IRAS positions are given in italics). VLA positions arc conservatively estimated to be accurate to 0.3" and IRAS 
positions to 20". (5): Heliocentric velocity (cz) from NED. (6): Luminosity dist ance (see Section 2). (7): FIR luminosity 
(see Section 3.3). (8): 1.4 GHz luminosity calculated from the NVSS catalog l lCondon et al.||1998l) . (9): FIR spectral 
luminosity (see Section 3,3). (10)): q,^ — log(LFiFt/Li.4) . 
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Table 3 (Full Version) 
Radio and optical observations and results 



IRAS VLA VLA S4.8 total 94.8 54.8000,5 0Mxe,„ logTj, VLBI Optical 

Name Array Morph (mjy) (mJy) I" ^"] [K] Epoch Class 

(1) (2) (.■?■) (4) (5) (f,-) (7) (H) (9-) (10) 





BnA— A 


CE 


5 66 


2 97 


00 






1* 


g 




BnA-A 




18.40 


2.85 


1.23 




2.63 










CE 


31.60 


2.95 


13.60 


1 97vn fifi 


2.75 


2 




lju u ir J. ^ u 


BnA 


E 


6.07 


3.13 


3.40 


Q 1 c; y-n Q1 


1.61 


\ 




nn'^4S-l-4'i'^ 1 




CE 


13 21 


2 95 


13.11 


1 1^1 V 41 


2 87 


2 




















2 




ni c,n'^_i_i 997 

y) ±0\.I0 1 


BnA 


CE 


6 23 


3 18 


5 25 




2 05 


1 






\ 


CE 


4 38 


3 23 


3 92 


1 13x0 74 


2 21 


1-2 






BnA 


CE 


8 23 


2 92 


1 85 


9 04 V 1 Ofi 


1.47 




L 




BnA 


CE 


11 16 


2 95 


7 37 


1 ^9v1 09 


2 28 


^ 








CE 


2.49 


3.39 


1.82 




2.45 


2 


L 






E 


12.20 


3.15 


0.93 




1.83 


2 




090S0-l-'^79f^ 




u 












1-2 




091 ^9-H 41 S 


BnA 


CE 


2.39 


3.82 


1.34 




1.52 


1-2 




n99ns -1-4744 




CE 


10 92 


3 08 


10 61 


1 7>4 y 1 09 


2 31 


2 




099'n'^-l-l Q99 


BnA 


u 












1-2 




^9'^4'^-l-9^'^'^ 

U ^ ^ ^ U vJ 


BnA 


CE 


14 69 


2 99 


3 61 


1 QO V 1 90 


1 74 


1-2 




09'^4fi-l-^41 9 




CE 


4.88 


3 15 


4 06 


9 1 fi V 1 9S 

ir. -LU A J- .irO 


1.71 


1-2 




09'^Q'^-l-'^4'^^ 






3.62 


3.34 


2.98 


5'^vO 4.0 


2.69 


2 




094^'=i-i-i 9"=;^ 


BnA 


CE 


20.60 


2.76 


14.98 


A 91 y 1 44 


1.93 


1-2 


g 


094^S_l_91 99 


BnA 




17 44 


2 61 


17 40 


98 V 1 


4 33 


2 




09'^nQ-l-l 94S 


BnA 


CE 


11.38 


3.05 


3.16 


9 1 1 V 1 ^9 


1.59 






^9'^1 1-1-1 9'^S 


BnA 


CE 


5 34 


3 07 


3 26 


1 Q1 V 1 '^1 


1 59 


1-2 




I ^ ^ tJ 1^ ^ U f 


BnA-A 


CE 


4 65 


3 35 


7 06 


fiO V 


2 88 










CE 


28 16 


2 68 


30 65 


9 1 vO fiQ 

ir. XU A U.Uc' 


2 86 


2* 


g 


n9'^79-^7^^9 




CE 


3 78 


3 30 


1 82 


R1 vO fi9 


2.10 


2 




O'il fi4-l-41 1 Q 




Q 


4Q4on no 


-0.76 


0.00 






1-2* 


L 




A 


CE 


4.79 


3.19 


1.02 


1 07x0 24 


2.14 


2 


H 


03266+4139 


A 


u 
















03406+3908 


A 


c 


0.40 


4.24 


0.41 


0.26x0.00 


3.04 


2 


S 


03449+7252 


A 


CE 


4.98 


3.22 


2.12 


0.81x0.67 


2.13 


2 




03514+1546 


BnA 


CE 


4.07 


3.27 


3.51 


1.56x0.93 


1.92 


1 




04002+0149 


BnA 


CE 


5.99 


3.10 


7.80 


4.68x1.78 


1.51 




L 


04007+2201 


BnA 


CE 


2.98 


3.39 


1.31 


1.02x0.80 


1.74 


1 


C 


04097+0525 


A 


U 














H 


04356+6738 


A 


CE 


6.07 


3.00 


5.61 


1.53x1.25 


2.01 


2 




04435+1822 


BnA-A 


C 


7.90 


2.91 


4.64 


1.88x0.05 


3.23 


1-2* 


H 


04520+0311 


BnA 


CE 


11.04 


2.93 


4.42 


3.42x1.25 


1.55 




C 


05054+1718 


BnA-A 


C 


9.45 


3.14 


7.43 


0.33x0.25 


3.49 


1* 


s 


05091+0508 


BnA 


cer 


13.06 


2.87 


0.43 


1.15x0.67 


1.29 




C 


05134+5811 


A 


CE 


1.62 


3.64 


0.57 


0.27x0.06 


3.09 


2 




05179+0845 


BnA 


CE 


3.68 


3.57 


1.30 


1.68x0.97 


1.44 




H 


05336+5407 


A 


CE 


5.48 


3.06 


3.56 


0.91x0.34 


2.60 


2 




05365+6921 


A 


C 


1.36 


3.97 


1.45 


0.06x0.00 


4.22 


2 




05405+0035 


BnA 


CE 


4.99 


3.34 


5.53 


1.42x0.87 


2.19 




H 


05414+5840 


A 


CE 


21.80 


3.01 


15.75 


0.83x0.39 


3.23 


2 




06052+8027 


A 


U 












2 




06140+8220 


A 


c 


57.69 


2.02 


58.72 


0.20x0.04 


5.41 


2* 




06239+7428 


A 


c 


1.69 


3.60 


0.00 






2* 




06538+4628 


A 


CE 


6.11 


3.24 


3.97 


0.57x0.40 


2.78 


2 




07062+2041 


BnA 


CE 


6.77 


3.15 


2.48 


0.97x0.81 


2.04 


1 


H 


07063+2043 


BnA-A 


CE 


7.07 


3.32 


4.55 


1.18x0.72 


2.27 


1 




07204+3332 


A 


CE 


7.77 


3.06 


1.68 


0.56x0.42 


2.39 


1-2 




F07258+3357 


A 


CE 


12.95 


3.23 


13.26 


0.56x0.35 


3.37 


1-2* 


H 


07258+3357 


A 


U 












1 




07329+1149 


BnA 


C 


10.62 


3.05 


10.23 


0.75x0.55 


2.93 


1 


H 
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Table 3 (Full Version) — Continued 



IRAS 


VLA 


VLA 


54.8 total 94.8 


54.8 comp 




logTb 


VLB! 


Optical 


Name 


Array 


Morph 


(mJy) 




(nijy) 


["x"] 


[K] 


Epoch 


Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


07336+3521 


A 


U 






- 


- 


- 


1 


- 


07566+2507 


BnA 


CE 


3.57 


3.27 


3.25 


1.65x1.07 


1.80 


1 


C 


08287+5246 


A 


CE 


4.58 


3.17 


1.21 


0.91x0.45 


2.01 


2 


- 


08339+6517 


A 


CE 


6.04 


3.08 


1.69 


0.62x0.48 


2.29 


2 


H 


08354+2555 


BnA-A 


C 


71.49 


2.61 


72.08 


0.44x0.30 


4.28 


1* 


- 


08561+0629 


BnA 


CE 


2.26 


3.39 


1.93 


2.01x1.32 


1.40 


- 


H 


09437+0317 


BnA 


U 






- 


- 


- 


- 


- 


10195+2149 


BnA 


U 






- 


- 


- 


1 


- 


11413+1103 


BnA 


CE 


8.02 


2.88 


1.13 


2.58x0.78 


1.29 


1 


H 


11547+2528 


BnA 


E 


8.68 


3.03 


9.68 


4.81x0.78 


1.95 


1 


- 


12099+2926 


BnA 


E 


11.89 


2.84 


13.26 


4.43x1.53 


1.83 


1 


- 


12159+3005 


BnA 


C 


14.16 


2.55 


14.24 


0.36x0.26 


3.72 


1 


S 


12208+0744 


BnA 


CE 


10.08 


2.85 


2.37 


1.37x0.56 


2.03 


2 


H 


13126+2452 


A 


C 


19.96 


3.03 


20.23 


0.20x0.14 


4.40 


1* 


- 


13183+3423 


A 


CE 


49.12 


2.65 


39.13 


0.58x0.29 


3.91 


1-2* 


C 


13188+0036 


BnA 


C 


7.64 


3.11 


7.13 


1.15x0.78 


2.44 


2 


C 


13238+3611 


A 


CE 


6.36 


3.02 


1.84 


0.75x0.52 


2.21 


1-2 


- 


13373+0105 


BnA 


U 






- 


- 


- 


- 


- 


13470+3530 


A 


c 


3.87 


3.42 


3.82 


0.27x0.07 


3.84 


1-2* 


- 


14003+3245 


A 


CE 


9.48 


2.98 


2.37 


0.77x0.51 


2.32 


1-2 


c 


14221+2450 


A 


C 


8.40 


2.92 


5.37 


0.35x0.29 


3.26 


- 


c 


14280+3126 


A 


U 






- 


- 


- 


1-2 


H 


15005+8343 


A 


CE 


5.28 


3.13 


1.30 


1.04x0.49 


1.95 


2 


- 


15107+0724 


BnA 


C 


31.69 


2.94 


32.34 


0.30x0.26 


4.16 


- 


c 


15276+1309 


BnA-A 


CE 


11.23 


3.06 


3.40 


0.64x0.51 


2.56 


1 


c 


15327+2340 


BnA 


CE 


183.13 


2.85 


97.51 


0.21x0.14 


5.06 


1* 


L 


15437+0234 


BnA 


CE 


6.69 


3.30 


6.51 


1.99x1.67 


1.83 


- 


s 


16030+2040 


BnA 


CE 


14.59 


2.78 


1.23 


0.61x0.31 


2.35 


1 


H 


16478+6303 


A 


E 


41.48 


2.20 


36.32 


2.48x1.43 


2.55 


2* 


- 


16577+5900 


A 


C 


15.24 


2.98 


7.80 


0.25x0.21 


3.71 


2* 


C 


17180+6039 


A 


U 






- 


- 


- 


2 


- 


17468+1320 


BnA-A 


c 


24.43 


2.57 


16.07 


0.17x0.13 


4.40 


1-2* 


- 


17530+3446 


A 


u 






- 


- 


- 


1-2 


- 


17548+2401 


BnA-A 


c 


16.19 


2.73 


12.08 


0.70x0.33 


3.26 


1* 


L 


18131+6820 


A 


CE 


6.37 


3.18 


6.34 


2.19x0.59 


2.23 


2 


C 


18145+2205 


BnA 


CE 


6.80 


3.13 


6.30 


2.67x0.60 


2.13 


1 


- 


18263+2242 


BnA-A 


CE 


3.73 


3.23 


1.15 


0.87x0.52 


1.94 


1 


C 


18425+6036 


A 


CE 


15.37 


3.01 


2.29 


0.13X0.00 


4.09 


2 


c 


18495+2334 


BnA 


CE 


6.43 


2.93 


6.95 


1.76x1.61 


1.93 


1 


c 


19000+4040 


A 


CE 


2.77 


3.53 


1.38 


1.25x0.44 


1.94 




H 


22025+4205 


A 


C 


15.72 


2.92 


13.01 


0.48x0.28 


3.53 


2 




22171+2908 


BnA 


C 


6.40 


3.17 


4.28 


1.33x0.55 


2.31 


1 




22388+3359 


A 


CE 


10.86 


3.04 


7.03 


0.51x0.42 


3.06 


1-2 


C 


23007+0836 


BnA-A 


C 


63.66 


2.75 


46.20 


0.33x0.25 


4.29 


2* 


s 


23106+0603 


BnA-A 


CE 


6.15 


3.00 


3.79 


0.45x0.22 


3.12 




H- 


23157+0618 


BnA-A 


C 


12.04 


2.94 


12.05 


1.02x0.52 


2.90 


2 


s 


23213+0923 


BnA 


CE 


5.34 


3.10 


2.27 


1.73x1.36 


1.52 




H 


23262+0314 


BnA 


CE 


7.99 


3.12 


1.52 


1.38x0.58 


1.82 




s 


23414+0014 


BnA-A 


CE 


10.86 


2.78 


3.77 


0.51x0.23 


3.05 


2 




23445+2911 


BnA 


U 












1 




23485+1952 


BnA 


c 


0.54 


4.40 


0.51 


0.86x0.40 


1.71 


1 


c 


23488+1949 


BnA 


CER 


23.88 


3.09 


3.30 


1.67x1.00 


1.84 


1* 


L 


23488+2018 


BnA-A 


C 


31.23 


2.87 


20.77 


1.17x0.73 


2.93 


1-2* 


C 


23591+2312 


BnA 


C 


14.09 


2.85 


2.30 


1.63x0.57 


1.93 


1 


c 


rUMNS — (1): IRAS Name 


(2): VLA 


array (s) 


used (3): 


Morphology class of VLA 4.8GHz 


image. 


C-Compact, 



Extended, CE-Compact plus Extended, CER-Compact plus extended with ring, U-Undetected. (4): Recovered 4.8GHz 
flux density in VLA image. (5). Log of ratio of 4.8GHz recovered VLA flux density to FIR flux density. (6): Total 
flux density of gaussian fitted to brightest feature in VLA imago {after removing VLBI component, hence zero indicates 
that brightest feature was dominated by the VLBI component). (7): Fitted gaussian major and minor axes to peak 
VLA feature (FWHM arcs ec). Note for IRAS15327+2340 (Arp 220) these are for the western nucleus as taken from the 
literature (sec Scction |4.5| l. (7): Fitted gaussian major and minor axes to peak VLA feature (FWHM arcsec). (8): Peak 
4.8GHz brightness temperature of gaussian component (see Section |4.5t . Note if formally unresolved along one axis a 
50mas FWHM was assumed along that axis). (9): VLBI Epoch(s) observed (asterisk indicates VLBI detection, see Table [4] 
for details). (10): Optical Class (sec Scction ftTet with H^HII, S^Scyfcrt, L^LINER and C^Composite. 



27 



- 001555+0250 0.66 - 



I 



01556+2507 3.23 




01579+5015 1.05 



_Q3 



- 02071+3857 0.36 - 



02152+1418 0.59 



02208+4744 1 .55 



93 



02345+2053 1 .06 







r02946+3412 


0.33 ; 

















02435+1253 2.42 




02395+3433 1 .76 



02509+1248 0.85 




3.2. — See caption of Figure 15^1-3.71 
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29 



05134+5811 0.48 




05179+0845 0.53 



<5) 



05336+5407 


1.51 











05365+6921 1.47 



05405+0035 2.40 



05414+5840 6.40 



I 



06538+4628 2.03 



07062+2041 1 .47 



07063+2043 0.70 



07204+3332 0.96 




r 



07329+1149 7.50 



(3 



07566+2507 1 .27 



3.4.— See caption of Figure 
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31 



14280+3126 



6.33 



& 



- 15005+8343 


0.46 







15107+0724 29.81 



15276+1309 


1.33 









15437+0234 2.17 



16030+2040 1.04 



140^1 



18131+6820 0.94 



18145+2205 1.98 



18263+2242 


0.38 






if 


4 






It 



- 18425+6036 2.19 - 



1 8495+2334 1 .48 




19000+4040 0.42 



3.6. — See caption of Figure 15^1-3.71 
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22025+4205 8.26 



22171+2908 2.40 



22388+3359 4.28 



06+0603 2.63 



23157+0618 3.87 



23213+0923 0.69 



23262+0314 0.91 - 



' 1 



23414+0014 2.62 



23485+1952 0.49 



23591+2312 1.14 



3.7.— See caption of Figure 
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